Dynamics of cholesteric structures in an electric field 
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Motivated by Lehmann-like rotation phenomena in cholesteric drops we study the transverse drift 
of two types of cholesteric fingers, which form rotating spirals in thin layers of cholesteric liquid 
crystal in an ac or dc electric field. We show that electrohydrodynamic effects induced by Carr- 
Helfrich charge separation or flexoelectric charge generation can describe the drift of cholesteric 
fingers. We argue that the observed Lehmann-like phenomena can be understood on the same basis. 
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One of the most interesting manifestations of macro- 
scopic chirality are cholesteric liquid crystals. In 
cholesterics, as in nematic liquid crystals, there is long- 
range orientational order of the elongated molecules 
along a local axis described by the director n. Whereas 
in nematics the elastic forces (or torques) tend to es- 
tablish a uniform orientation of n, the chiral molecules 
in cholesterics lead in equilibrium to a helical arrange- 
ment with n perpendicular to the helix axis. Choosing 
the axis along z, the structure with pitch po is given by 
n = (cos tp(z), sin <p(z),0), <p = qoz, qa = 2w/p . 

The helical symmetry leads to interesting dynamical 
effects. Particularly intriguing is the Lehmann rotation 
of the director structure in a cholesteric droplet heated 
from below, described in 1900 Q, SEI- ^ nas not been 
observed again, which is usually attributed to the in- 
fluence of surface anchoring. However, the electric ana- 
log where the temperature gradient is replaced by a dc 
electric field E ( "electromechanical effect" ) has been ob- 
served ,4j. Traditionally, the explanation is based on 
phenomenological hydrodynamic considerations, which 
by symmetry allow for an additional dissipative dynamic 
coupling between the director and electric field ("elec- 
tromechanical coupling") 0,0]. Strictly speaking, here 
"hydrodynamic" means that spatial modulations must be 
slow on the scale of the pitch which is not the case in Q . 
In fact, in the experiments the director structure inside 
the droplet was not a homogeneous helix, but included 
substantial spay-bend distortions and even defects. An- 
other weakness of the approach is that no underlying 
mechanism has been identified that would, at least in 
principle, allow to determine the coefficients involved. 

We have investigate for the first time a driving mecha- 
nism for chirality-related dynamical phenomena, involv- 
ing well-established electro-hydrodynamic (EHD) effects. 
Actually, director rotation is not the best choice to study 
such phenomena, since special precautions are needed to 
avoid surface anchoring Another effect is related to 
the intrinsic length scale defined by helical symmetry. In 
a confined geometry this can lead to the spontaneous for- 
mation of spatial structures ("cholesteric fingers") lack- 
ing certain reflection symmetries, which under nonequi- 



librium conditions results in drift of the structures. This 
is not affected by surface anchoring. We will show that 
this drift can be explained by EHD effects and, at the 
end, return to the rotation phenomena. 

Cholesterics placed between two plane parallel elec- 
trodes (separation d) , providing strong homeotropic (per- 
pendicular to the electrodes) anchoring of n, can experi- 
ence unwinding of the helix due to the orienting effect of 
the electric field and of the boundaries. The unwinding 
transition, which is typically discontinuous, occurs when 
a combination of the confinement ratio C = d/po and 
the electric field strength (or applied voltage U) reaches 
a critical value [5j. Thus, there exists a line in the (U, C) 
plane where the two phases coexist. Near this line one 
finds the cholesteric fingers (CFs): elongated structures 
that are localized or arrangedperiodically. At least four 
types of CFs were observed [fjj. The director configura- 
tion of CF of the first type (CF1) is invariant with respect 
to a 7r- rotation about the finger axis , whereas CF2 has 
a mirror symmetry with respect to the midplane of the 
cell see Figs. E] ED (top). The difference in structure 

of CF1 and CF2 manifests itself in the dynamics. In a dc 
field both fingers are observed to drift perpendicular to 
their axis [9j, Il0| , whereas in an ac field only CF2 drifts 
0, 0, 0| . Since CFs also grow along their axis one ob- 
serves the formation of CF spirals. Most measurements 
of the transverse drift of CFs are based on the analysis 
of spiral dynamics. 

As the motor of the CF1 drift only the electromechan- 
ical coupling has been proposed [lC| . The magnitude of 
the drift velocity estimated with the coupling coefficient 
taken from droplet rotation-type experiments turns out 
to be at least one order higher than observed [T^j. To 
explain the drift of CF2 in an ac field, several models 
were proposed 0, including electromechanical coupling. 
However, these models fail to describe the recently ob- 
served falloff of the drift velocity V± when the frequency 
/ = u>/2tt of the applied electric field approaches the in- 
verse charge relaxation time T q [l2T ]. In the experiments 
the conductivity, and thereby r q , was varied by using 
different concentrations of ionic dopant. 

The mechanism we propose is based on flow induced ei- 
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ther by charge separation through anisotropic conductiv- 
ity, i.e., the Carr-Helfrich effect, (CF2 under ac driving) 
or by flexoelectric charge generation (CF1 under dc driv- 
ing). We use the standard set of nematodynamic equa- 
tions for the director n, the velocity v (Navier-Stokes 
equation) in the Stokes approximation (neglect of iner- 
tial terms), and the electric field E |2j, |3j. As usual 
the electric properties of the material are described by 
a dielectric permittivity tensor ey = ex#y + £ a ninj, 
(ex, £|| = £_l + e a are the permittivities perpendicu- 
lar and parallel to the director, respectively), the anal- 
ogous conductivity tensor cry = a±6ij + cr a mnj, and 
a flexopolarisation P? 1 = ein(V ■ n) + 63(71 • V)n. 
We introduce dimensionless variables t — r^t, r = df, 
E = (U /d)E, P-f 1 = (e Q U /d)P fl , and charge density 
p e i = (e U /d 2 )p e i. Here r d = 71 d 2 / K 33 is the relevant 
director relaxation time (71 = rotational viscosity, -K33 = 
bend elastic constant), Uo = \J K 33 / eo characterizes the 
typical voltage (Uo « IV for the materials used). 

We choose the z axis perpendicular to the bounding 
electrodes and allow for a drift of the structure with ve- 
locity V± in the x direction, transverse to the finger's long 
axis by replacing dt — > dt — V±d x . Then the equations 
can be written as (tildes are omitted, V± = V±x) 

[d t + (v - V±) ■ V + iztA- Ox]n = -l ± h r , (1) 
p ti - T] % 3 + h v k n k , t = p el Ei - V±Sjij , (2) 

^[d t + {v-v^)-v] Pel + p el = 

Td 

V- [-e^ H (n-E)n + P fl ]. (3) 

supplemented by the incompressibility condition V • v = 
0, the Poisson equation p e i = V- [ei_E+e a (n-E)n+Pf l ], 
which was already used in the charge conservation equa- 
tion ©, the electrostatic condition Vx£ = 0, and the di- 
rector normalization n 2 = 1. The notation = df/dxi 
is used throughout. The generation of space charges is 
characterized by the Hclfrich parameter = cr a /<7j_ — 
e a /e x in Eq. ©. 

The director equation couples to the flow field 
through the local fluid rotation = (V X v)/2 and the 
hydrodynamic strain tensor Ay = (djj + with 
72 = 72/71 and the projection tensor <5y = <5y — mrij. 
Coupling to the elastic and electric torques is through 
h\ = SF/Srii with the free energy density F = |fci(V • 
n) 2 + \k 2 (2irC + n • (V x n)) 2 + ±fc 3 (n x (V x n)) 2 - 
±e a (n • E) 2 - P fl ■ E. Here k t = Ku/K 33 with the elastic 
constants Ka. 

In Eq. J2J) the Stokes approximation is justified since 
the processes of interest are controlled by Td ~ 1 s and 
the charge relaxation time r q = eQe±/a± ~ 10~ 3 s, 
which are much larger than the viscous relaxation time 
Tu = p m d 2 /ji ~ 1CP 6 s (p m = mass density). The elas- 
tic part of the stress tensor (Ericksen tensor) has been 
eliminated and the pressure redefined: p = po + F |l3| . 
The viscous stress tensor is Xy = d x riinjn}-n m A\~ m + 



a' 2 riiNj + a' 3 rijNi + a^A^ + a' 5 niii k A kj + a' 6 njn k A ki , 
where o! i = ctj/71 with the Leslie viscosity coefficients 
Qij. Moreover N = (dt + v ■ V)n — ft x n and h k = 
N k +^' 2 njAj k . The bulk force p e iE is the Coulomb force 
and = a 2 nid x rij + a 3 rijd x ni. 

We solve Eqs.©-© in a perturbative way: n = n + 
m + ...,E = E + Ei + .. ., v =wi + ..., V± = V±i + ..., 
Pel = Peii + • • •• At lowest order the electric charge is 
neglected. Therefore no bulk force arises in Eq.© and 
thus i>o = 0. no and Eq are obtained from 

S^h^O, V-[e x E + e a (n -E )n ] = 0. (4) 

We write Eo — Eq(z — V0o) cos(uiTdt) with Eq the ap- 
plied electric field and <po the induced potential. Previ- 
ous results show that the director configurations obtained 
from Eqs.(@J and the stability diagrams are in good agree- 
ment with the experiments 0, la, 0|. We have solved 
Eqs.QJ in two dimensions (infinite extension of the CF 
in the y direction) by a relaxation method. 
At first order Eq.© gives 

cos(ujT d t) + LUT q sm(wTdt) t 
Peii — ^— — n o -&oe±4-ff 

xV- {[n • (z- V0o)]«o} 

+ e/ ;V-{no(V-n ) + — (n -V)n }, (5) 
ei 

where e/z = e\d/ (eoUo)- The dc case is covered by set- 
ting u> = 0. Next one solves the Navier-Stokes Eq.© (to- 
gether with the incompressibility condition) at first order. 
Taking the curl eliminates the pressure and leads to an 
inhomogeneous ODE for V\ with the coefficients depend- 
ing on no- The form of the inhomogeneities [right-hand 
side of J5J|] give solution of the form 

e £ E 2 

vi = ^" 2 e/ / Si + e n Eoh + Vli/3, ( 6 ) 

where -Ee// = E /V2 (Eo) for ac (dc) driving and the 
second term appears only in the dc case (then also u) = 0). 
For the ac case a time average of Eq.© was taken. The 
functions fi depend only on no and the a[ (f% depends 
also on the ratio of the flexocoefficients 63/ e±). We have 
determined the functions fi [from the linearized Eq.©] 
by a Galerkin expansion with trigonometric functions (for 
CF1) and Hermite polynomials (for CF2) as trial func- 
tions. 

Finally, the director equation at first order gives 

S^hl + S^hl = (vi-V-V X id x )n 

+^'26^-^710 - Oi x n , (7) 

where h^, 5i are linear in n\. The homogeneous prob- 
lem <5p_ hi +5i_ ± hQ = is solved by the translation mode 
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FIG. 1: Director profile (top) and induced velocity profile 
(stream lines) (bottom) of CF2. U = 1.9V, lot„ = 1, 5CB 
material parameters (l4| |. 



16 



12 -- 



> 



4 - 




10 



* «> 


i— 


* n 


_ U X 

o T JT''>-\ 





















10" 



10"' 
ore 



10 



10' 



FIG. 2: Drift velocity of CF2 versus reduced frequency: 
points are experimental data from |T3l , lines are calculations 
from JSJ: for a\\/(r± — 1.3 (solid), a^/a± — 1.4 (dashed), and 
CT||/crx = 1.5 (dash-dotted). C = 1.77, d = 31 /im, f/o = 1.9 
V, 5CB material parameters p^| . 



dno/dx. The solvability condition for the inhomoge- 
neous problem [with V\ substituted from ©] fixes the 
drift velocity V±±. In physical units we obtain 



T/ d e e±£_ H r _> 

7i 1 + w 2 r 2 e// Jo 



/a 



e\E Q I 2 
7i + ^3 



(8) 



where J = (n 0:X -rio^), 7i = (n , x -Si), 7 2 = (n ,x -92), 
73 = (no,i • 93). The last term appears only in the dc 
case. The functions g±, g 2 , gz are easily expressed in 
terms of the fi, f 2 , /3, respectively. The scalar product 
is defined by (a ■ b) = J J (a ■ b)dxdz. 

We first discuss our results on CF2 drift in an ac field. 
Since in the experiments the arms of CF2 spirals are well 
separated [l2^ , we have in the computations chosen a box 
width large compared to the width of the director struc- 
ture (isolated finger). The director profile n together 
with the stream lines of the flow from 10 as computed 
for the parameters of the material 5CB (4'-n-pentyl-4- 
cyanobiphenyl) |l4j used in the experiments are shown 
in Fig. ^ 111 Fig. [21 our results for V±i, which is of 
order of v\ , versus reduced frequency u>r q are shown, to- 
gether with the experimental results (the different sym- 
bols relate to different impurity concentrations between 
2 ■ 10~ 5 and 0.05 wt%) .12]. The authors of pjj have 
scaled the frequency down by a factor of about 1.75 in 
order to account for the fact that the charge relaxation 
time T q was measured in the isotropic phase at temper- 
ature T = 40 °C (the experiments were done at T = 30 
°C). The three curves correspond to a\\/a±_ = 1.3, 1.4 and 
1.5, which agrees well with measurements in pure 5CB 
where cr\\/a± = 1.44 0. <J\\/a± is the only parameter 
not given in [l^ . 

Another experimental result born out by the model is 
the approximate linear dependence of the drift velocity 
of CF2 as a function of applied field for samples with 
different thicknesses but fixed confinement ratio C \\'M. 



The CF2 appear only around an electric field such that 
(E d) 2 = 7r 2 /(e e a )(AC 2 K 2 2 - K 2 3 )/K 33 @. Using this 
to eliminate d from © one obtains the linear dependence 
of the drift velocity on Eq . 

We now turn to CF1 in a dc field. From the symme- 
try of the director profile n one now has I\ = 0, so 
one is left with the last term in Eq.® relating to flex- 
oelectric charge generation. Unfortunately one now has 
to cope with various uncertainties. First, the approxi- 
mation of isolated fingers is not valid, since in the ex- 
perimental spirals neighboring fingers are packed closely 
JloL Hr]| . Second, in the dc case one has to expect 
screening of the electric field by Debye layers, as evi- 
denced in 0| and also suggested by the voltage offset 
~ 2V in the current-voltage curve presented in [lQj [both 
using the material MBBA (4-methoxy-benzylidene-4'-n- 
butylaniline)] . Thus we have solved Eqs.Q on a width 




FIG. 3: Director profile (top) and induced velocity profile 
(stream lines) (bottom) of CF1. U = 0.2V, L = Lf = 1.95, 
MBBA material parameters |T?|. 
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L in the x direction with periodic boundary conditions. 
Vj_i turns out to be sensitive to L: it decreases with 
decreasing L and even changes sign at L sa 1.6. In the 
absence of experimental data on the finger width we have 
- on one hand - minimized the free energy density with 
respect to L for given values of the electric field leading 
to the "optimal" box width Lp(U). Typical values of 
Lp (in units of d) are 2 2.5 for a voltage 0.2 3V. In 
Fig. an example of the calculated director profile and 
the stream lines of the flow from © are shown. We have 
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voltage (V) 

FIG. 4: Drift velocity of CF1 versus voltage: points are 
experimental data from [To|| shifted by 2V; lines are calcula- 
tions from IjHJl. Solid line: L — Lf(U); dashed: L — 1.9; 
dot-dashed: L — 2.1. Flexocoefficients ei = —1.05 • 1CF 11 
C/m, e 3 = -1.25 • 10" 11 C/m. C = 1.14, d = 12 /im, MBBA 
material parameters |l7|. 

also made calculations for some fixed values of L. Our re- 
sults for V±i as a function of U, corrected by a screening 
of 2V, are given in Fig. 0] together with the experimental 
data |l(J . Typical values of the flexocoefficients for pure 
MBBA were chosen [li |. 

In conclusion, we have developed an EHD model for 
the drift of CF1 and CF2. For CF1 in a dc electric 
field, flow induced by flexoelectric charge generation can 
describe the drift, but no quantitative conclusion could 
be drawn. For CF2 a Carr-Helfrich-like mechanism de- 
scribes quantitatively the drift in an ac electric field. 
Moreover, our preliminary studies show that the EHD 
model with flexoelectric charge generation (as in CF1) 
can describe the rotation of cholesteric droplets in a dc 
electric field Q, which has hitherto been interpreted in 
terms of phenomenological electromechanical coupling. 
Furthermore, it seems very likely that analogous thermo- 
hydrodynamic effects, which lead to very efficient convec- 
tion phenomena in liquid crystals , can also describe 
the original Lehmann rotation [1|. This then suggests 
that so far there is no clear experimental manifestation 
of the (unspecified) phenomenological electro- or thermo- 
mechanical coupling. The EHD model provides a very 



general mechanism for forces and motion that can be 
applied to other director structures, like other types of 
cholesteric fingers (CF3 and CF4) and defects in an elec- 
tric field. 
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